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^ ■ Abstract 

^ ; 

' Symmetry-projected Hartree-Fock-Bogoliubov (HFB) equations are derived using 

the variational ansatz for the generalized one-body density-matrix in the Valatin 
form. It is shown that the projected-energy functional can be completely expressed 
J> I in terms of the HFB density-matrix and the pairing-tensor. The variation of this 

^2 ' projected-energy is shown to result in HFB equations with modified expressions 

. for the pairing-potential A and the Hartree-Fock field F. The expressions for these 

I quantities are explicitly derived for the case of particle number-projection. The 

numerical applicability of this projection method is studied in an exactly soluble 



o 

. model of a deformed single-j shell. 



^ . 1 Introduction 

> 

^ . The concept of spontaneous symmetry breaking in mean-field theories has 

^ ! played a central role in our understanding of many-body problems [1]. The 

exact solution of a quantum mechanical many-body problem is impossible 
except for a few-body system and it is therefore necessary to adopt approxi- 
mate methods. The most popular of these methods has been mean-field theory 
based on effective forces. The basic idea of the mean-field approach is to ap- 
proximate the "unknown" many-body wavefunction by a Slater-determinant 
of single-particle or quasiparticle configurations. In this way the many-body 
problem is reduced to an effective one-body problem. In other words, the 
system of interacting particles is transformed to a system of non-interacting 
particles or quasiparticles. The effective one-body potential is obtained from 
the many-body Hamiltonian by using the Hartree-Fock (HP) or Hartree-Fock- 
Bogoliubov (HFB) approaches [2]. The HF method has first been introduced to 
solve the atomic many-body problem and it was shown that the many-electron 
problem can be reduced to one-electron problem with an effective potential 
which is determined by the interactions among all the electrons. This method 
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has been successfully applied to the nuclear many-body problem, although at 
first sight it appears conceptually contradictory to use a mean-field concept 
to a strongly interacting system. It was established later [3] that the Pauli- 
exclusion principle favours the description of the nuclear many-body problem 
in terms of the mean-field concept based on effective interactions. 

The HF approximation describes the long-range part, referred to as the particle- 
hole (p — h) channel, of the effective interaction and, indeed, is quite appropri- 
ate for the atomic many-body problem. However, for the nuclear many-body 
problem it can only be used in closed shell configurations. For open shells, 
strong short-range pairing-correlations in the particle-particle (p — p) channel 
need to be considered. In oder to treat the p — p channel, the Bardeen-Cooper- 
Schrieffer (BCS) approximation [4] , originally introduced to explain supercon- 
ductivity in metalls, has been applied to the nuclear problem with remarkable 
success [5,6]. For a complete description, both p—p andp — h channels need to 
be considered in a self-consistent manner. This is achieved in the generalised 
Hartree-Fock-Bogloliubov formalism [7]. 

It is an essential feature of mean-field theories that they can violate sym- 
metries, i.e. the approximate product wavefunctions does not obey the same 
symmetries as the underlying two-body Hamiltonian. The broken symmetries, 
for instance, include translational symmetry, rotational symmetry, and the 
gauge symmetries connected with the particle numbers. Since the exact so- 
lution of the many-body Schrodinger equation obeys these symmetries being 
eigen-function of a complete set of generators of the corresponding symmetry 
group, the violation of symmetries in the first place is a very undesirable fea- 
ture. However, the concept of symmetry violation allows to take into account 
correlations in a very simplified form by employing product states which are 
connected with phase-transitions and provide a relatively easy understanding 
of very complicated phenomena. Nevertheless, in systems with finite particle 
number, such as atomic nuclei these phase-transitions are never sharp and are 
smeared out by fluctuations. It is one of the disadvantages of the mean-field 
approximation that it leads to sharp phase-transitions even in finite systems, 
where the exact solution which takes into account fiuctuations beyond the 
mean-field shows only a smooth transition between the two phases. 

In order to have a better description of finite systems in the vicinity of phase- 
transitions, it is necessary to include fluctuations. One very powerfull way to 
consider such fluctuations is through the restoration of the broken symmetries 
by using the projection methods [8]. 

There exists a vast amount of literature on such projection methods at zero 
temperature (see, for instance, Ref. [2]) and at flnite temperatures (see, for 
instance, Ref. [9]). Since mean- field theories are variational theories, one can 
carry out the projection before or after the variation. It has turned out that 
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variation after projection (VAP) [10] is the appropriate tool that fulfills the 
variational principle and provides a self-consistent description of fluctuations 
going beyond mean-field. 

Although, the method of variation after projection has been known since 
more than thirty-years, the numerical solution of the corresponding varia- 
tional equations is relatively complicated. Therefore, a fully self-consistent 
variation after exact projection has been carried out so far only in a limited 
number of cases. Most of such calculations have been restricted to the varia- 
tion with respect to a few order parameters. Applications with variation after 
exact projection have been restricted to light nuclei [11] or to the case of vio- 
lation of particle number in BCS-theory (see, for instance, FBCS in Ref. [12]). 
Exact projection within full HFB-theory is possible by a search of the mini- 
mum in the projected-energy surface by gradient methods [13]. However these 
methods are numerically very complicated. So far they have been applied only 
to the case of number-projection in restricted spaces and for separable-forces 
[14]. In most of the applications, where number-projection is necessary, the 
approximate method of Lipkin and Nogami is used, which is, in fact, only an 
approximation and violates the variational principle [15]. 

The solution of the unprojected variational problem can be represented as a 
non-linear diagonalization problem (HF- or HFB-equations). This allows an 
iterative solution by the application of very fast numerical techniques. So far 
the projected variational problem has not been written as a diagonalization 
problem. In this work we derive for the first time a set of symmetry-projected 
HFB equations. They have same form as in the unprojected case and the only 
thing which differs is that in the case of projection the expressions for the 
Hartree-Fock potential and the pairing-field are more involved. This allows 
to modify the existing codes of unprojected solutions of the HFB-equations 
without much difficulity. 

The manuscript is organised as follows: In section 2, some basic HFB relations 
are presented. The projected HFB energy is expressed in terms of the norm- 
and the Hamiltonian-overlaps in section 3 and it is shown in section 4 that 
these overlaps are entirely expressible in terms of the density-matrix p and the 
pairing-tensor k. In section 5, it is shown that the variation of an arbitrary real 
energy functional, which can be completely written in terms of the density- 
matrix and the pairing-tensor, results in HFB equations of the conventional 
form. The exphcit expressions of the projected HFB fields are derived for the 
case of particle-number projection in section 6. In section 7, it is shown that 
the expression for the pairing-field reduces to the familiar form of FBCS in the 
canonical representation. A model study is presented in section 8 and finally 
the present work is summarised in section 9. 
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2 Basic HFB-Formulae 



In second quantization, the many-body Hamiltonian of a fermion system is 
usually expressed in terms of a set of annihilation and creation operators 
(c, c^) — {ci, . . . , Cj^] Ci, . . . , Cj^): 



nin2 



1 — t t 

^ nin2n3n4 



where the anti-symmetrized two-body matrix-element is defined by 



nin2n3n4 



(2) 



The operators annihilate the bare fermion vacuum 

Cnl-) = 



(3) 



for all n and provide a basis in the present analysis. 

In the next step, a set of quasiparticle operators (a, a'^) = {a^, . . . , aj^; a\, . . . , Q;j\^) 
is introduced, which are connected to the original particle operators by a linear 
transformation 



n 

n 

This can be written in matrix form as 





' t/t yt ^ 




= w* 
















V) 







(4) 
(5) 



(6) 



Quasiparticle operators have to satisfy the same fermion commutation rela- 
tions as the original operators. Therefore, the transformation matrix is re- 
quired to be unitary 

which leads to following relations among the coefficients Unk and 
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U^U + V^V = I, UU^ + V*V^ = /, (8) 
U^V + V^U = 0, C/0 + V*U^ = 0. (9) 

The quasiparticle operators annihilate the quasiparticle vacuum |$), defined 

by 

= 0, (10) 

for all k. In mean-field theory, it represents an approximation to the ground- 
state of the system and turns out to be a generalized Slater-determinant [2]. 
Since the quasiparticle transformation mixes creation and annihilation oper- 
ators, 1$) does not correspond to a wavefunction with good particle- number. 
We therefore have two types of densities, the normal density p, and the pairing- 
tensor defined as 

Pnn' = ($|4c„|$), Knn> = ($|c„,C„|$). (11) 

These can be expressed in terms of the HFB coefficients as 

p = V*V^, K = V*U^ = -UVl (12) 

It can be immediately seen that p is Hermitian and k is antisymmetric, i.e. 

p^ — p, — —K. (13) 

Using the unitarity relations (7), it can be also shown that p and k satisfy 

P — P^ — —KK*, pK = Kp* . (14) 

With these basic relations for p and /t, the generalized density-matrix in the 
Valatin form [16] is defined as 



=K =wr" wt, (15) 







( p K ) 




(oo\ 


r \ 






lOll 


\ -K* -a* i 




[-K*l-p* 







which satisfies, 

= n. (16) 

It turns out that there is a one-to-one correspondence between the HFB- 
wavefunction |$) and the generalized density-matrix TZ, and that the quasi- 
particle transformation W diagoualizes this matrix. Using Thouless theorem 
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and restricting the following considerations to the case of even number-parity, 
we can express |$) in terms of the operators (c, c^) as 

\^)^Me^\-), (17) 

with 

Z^Y. Znn'cicl^ (18) 
n<n' 

where the matrix Z is given by 

Z = V*U*-^ = -pK*-^. (19) 

jV' is a normalization factor 

M= -\/|dett/| = det(l - Z*Z)-3. (20) 

3 The Projected-Energy 



In this section, we express the projected-energy in terms of the Hamiltonian- 
and the norm-overlaps. Restricting ourselves in this work to abelian symmetry 
groups, the projection-operator is defined as 

= jdgd\g)R{g), (21) 



where the integral runs over all elements g of the symmetry group. In the case 
of number-projection g — (f) is the gauge-angle, I — N is the particle-number 
and 

R{(f)) = e*^*, ci^(0) = -le-^-^^. (22) 

ATT 



For one-dimensonal angular- momentum projection of an axially-symmetric 
nucleus, g = j3 is the Euler-angle, / is the angular momentum and 

m^e'^'^, d\p)^^^dUP). (23) 



In the case of i^- projection in triaxial nuclei, i.e. projection onto a good angu- 
lar mometum component along the 2;-axis g = j is the azimuth angle around 
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the 2;-axis, I — K and 



R{l)^e^^'% rf^(^) ^ i_e-7i^. (24) 



Assuming that the Hamiltonian H commutes with the symmetry operator 
R{g)i the projected-energy is given by 



(*|P^|*) jdgd^{g)(m{g)\^) ' 



In order to simphfy the notation, we introduce the rotated Slater-determinants, 
l^) = 7^^^' ^ith |0) = |$) and (0|y) = l, (26) 



and the coefficients 

x{g) ^ d\g) U\R{g)\^) , (27) 



and 



These have the property 

x{-g)^x*{g), y{-g)^y*{g). (29) 



Here {—g) labels the inverse rotation R{—g) = R\g)- For the projected energy 
we then find 



The overlap functions (^^\R{g)\^J and {0\H\g) are evaluated in the following 
section. 
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4 The Overlap-Integrals 



It will be shown in this section that the overlaps for the norm {^\R{g)\^) 
and the Hamiltonian {0\H\g) can be completely expressed in terms of p and k 
defined in Eq. (11). For this purpose, we use the generalized Wick theorem as 
derived by Onishi (for details sec Rcf. [2]), which allows us to express overlaps 
in terms the HFB coefficients {U, V) and {Ug, Vg) corresponding to the Slater 
determinants |0) = |$) and \g) ~ R{g)\^), respectively The HFB-coefhcients 
{Ug, Vg) of the wavefunction \g) are given by 

U = DU V = D*V, (31) 



where the matrix Dg is the representation of the rotation operator R{g) in the 
single-particle space characterized by the basis states \n) — al\—) 

{Dg)^^, = {n\Rig)\n'). (32) 



This means 



Wg^VgW, with Vg^ 



^DgO^ 

d; 



(33) 



In order to derive the norm-overlap, we start from the representation (17) for 
the HFB-wave-functions and find 

($1^(^)1$) = I detC/| (-|e^^e^«|-) with Zg^V*U*-K (34) 



We introduce the function 



G,(A) = (-|e^'e^^«|-), (35) 



which obeys the differential equation 



_d „ ... 1 
dX 



Gg{X)^-Gg{X)Tr:\n{l-XZ*Z^). (36) 



Integrating from A = to A = 1 we find 

^.(A) = e^'^i°(^-^^'^«) = ±^det(l-AZ*Z^). (37) 
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The phase of the overlap integral remains open in this derivation, because 
the logarithmc is a multivalued function determined only up to multiples of 
2ni. Therefore, we are left with a phase problem. In principle, this phase is 
well defined. Several methods have been proposed in the hterature [17,18] to 
determine the phase in an arbitrary case. However, they appear complicated 
and are connected with considerable numerical effort. 

Using the expression for Gg{l) in (34), we obtain the norm-overlap 



($1^(^)1$) = I dett/|GA(l) = ±. 



det(l - Z*Zg) 
\ det(l - Z*Z) 



±JdetDjdet{Ng){38) 



with the matrix Ng defined as 



Ng = U'^U* + V'^V; = U^D*U* + V^DV*, 



(39) 



For the Hamiltonian-overlap functions, {Q\H\g), we again use the generalized 
Wick theorem of Onishi, which allows us to express matrix elements of arbi- 
trary operators between the brackets (0| and \g) for Slater-determinants \g) 
(normalized by {Q\g) = 1) in terms of the generahzed densities. The Hamito- 
nian overlap has three parts 



{0\H\g) = H{g) = H,,{g) + H„,{g) + H,,{g), 



E 



n\n2 
nin2 

4 n 



2 ^ Vn2n3n4^0|4iC„3l5)(0|cT^C„j5 



+ 7 E Vn2n3n4(0|4i42l^)\0|Cn4C„3|^ 
n\n2n3nA 



(40) 
(41) 

(42) 

(43) 



with 



Hsp{g) = E e„i„2Pn2ni(^) = 1^ {epig)) > 

n\n2 



Hph{g) 



nin2n3n4 



mn2n3n4Pn3n,(g)Pn4n2(g) = 2^ (^(g)p(g)) , 



nin2n3n4 nin2 



'^nin2ig)'^n3nA{g)j 



nin2n3n4 
1 



Tr(A(^)^*(^)) = -^Tr(A*(^)«:(^)), 



2 ' 2 
where we have introduced the fields 



(44) 
(45) 
(46) 
(47) 
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nina (9), (48) 

n2n4 

^711^2(5') ~o E '^nin2n3n4 '^713714(5')' (49) 

^71102(5') ^2 E '*n3n4(5')'^n3n4nin2' (^0) 
n3n4 

^^1712(5') = 2 '^nin2n3n4'^n3n4.i9)- (51) 

n3n4 

The generalized densities are given by 



Pnn'ig) = {0\clcjg) = V^)^^, , (52) 

^nn'ig) = (0|C„,CJ^) = {V;N;'U^)^^^ , (53) 

-<,,(^7) = {0\ci4,\g) = (f/;iV;V^)^^, , (54) 

-<n'{9) - (0|c„,4 |y) = . (55) 

The matrices K{g) and K{g) are anti-symmetric and the matrix a{g) is related 
to p{g) through 

aig)^p^{g)-l. (56) 



So far all the expressions still depend on the HFB coefficients U and V. Using 
the definitions for the densities p and k, we now rewrite expressions for the 
overlaps entirely in terms of these quantities as 



($|i?(c/)|$) = ± Jdet DgJdet{U^D*U* + V^D;V*] 



±Jdet Dg^det{V-\VUWluV^ + VVWJVV^){V^)-^} 



= ±ydet^(detp)"^/ydet (^pDgP - KD*gK*), (57) 

dot D„ 



wdetpdet Cg 
where we have introduced the rotated-densities 



(58) 



Ps-D^PDI P-g^DlpD^, (59) 
Kg^DgKDj, K.g^DlKD;, (60) 
Zg = D^ZD^, Z_g = DlZD; (61) 
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and 



Cg ^ = PPg- f^f^l- 



(62) 



The transition densities appearing in the Hamiltonian-overlaps can also be 
rewritten in terms of the HFB densities as 



It is now clear from the above expressions that the norm- and the Hamiltonian- 
overlaps are completely expressible in terms of densities p and n. In the next 
section, it will be shown that the variation of an arbitrary real energy func- 
tional which is expressible in terms of these densities results in the HFB equa- 
tions. 



5 Variational Equations 

We have now expressed the projected- energy as a functional of the HFB- 
densities p and k. So far we have used the fact that p in hermitian and that k 
is skew-symmetric, but we have not used the relations IZ^ = TZ {p — p^ = —kk* 
and pK — Kp*). Without these equations we have the independent variables 
^G{Pn'n),^'>TT'{Pn'n), Pnn, 3?e(/t„/„), and $5m(K„/„) for n < n' and we have to 
vary with respect to them under the constraint 7^^ = TZ. As shown in stan- 
dard text books, this variation can be replaced by a variation with respect 
to the independent variables Pn'n^ Pn'n^ Pnm i^n'm k*/^ for n < n'under the 
constraint, TZ^ — TZ. Introducing the Lagrangian multipliers A„„/ the varia- 
tional ansatz is 



As we have shown in the earlier section that the projected-energy depends 
only on p and k, the variation of the energy is expressed as 



p{a) = PgCgP, 

t^{9)=PgCgK, 

'^*{9)=f^*gCgP, 

a*{g) = K*gCgK. 



(63) 
(64) 
(65) 
(66) 




(67) 




+ 



dE 



^P*n'n + 




(68) 
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and introducing the quantities 
dE 



for n < n', A„„, 



dE 



for n < n'. 



Since the functional E is real, we find 

h* 



A* , 



dE 



It can be easily shown that 



and obtain 



dE ^ dE 



1 r as dE ^ ^ 

^ Vnn' ^Pn'n ^Pn'n 

— 21^ ^nn'^Pn'n + Km'^Pn'n ~ ^nn'^'^n'n ~ ^nn'^'^n'n^ \ ' 
V nn' ) 



1 



{Tr (Mp) + Tr {h*5p*) - Tr {A*Sk) - Tr {A5k*)} , 



It. 

2 



h A \ I Sp Sk 



-A* -h* j \ -Sk,* -5p 
and introducing the matrix 

/ 



h A 
-A* -h* 



we have 



SE - -Tr (nsn) 



Including the constraint leads to the variational ansatz 

5 {e{p, K)-Tr (A(7^^ -Tzj)}^ ^Tr {{H-ATZ -TZA + A)5TZ} 



12 



Since 5TZ is an arbitrary variation, we find 

n-ATZ - 7^A + A = 0. (79) 

Using (TZ^ — TV) , the above equation can be written as 

[n, n\ = 0, (80) 



which is solved by the HFB-equations 





(81) 



The HFB-equations have been obtained in (81) with the only assumption that 
the energy functional is expressible in terms of p and k. As we have shown 
that the projected- energy is also completely expressible in terms of these 
quantities. Therefore, in the case of variation after projection we find HFB- 
equations of the same structure. We only have to consider a projected-energy 
functional 



E' = \^i^pT^ - J dgy{g) m\g) , (82) 



which yields different expressions for the fields h and A: 



In the following section, these derivatives are obtained for the case of particle- 
number projection. The expressions of these derivatives in a general case be- 
come quite complicated and will not be presented here. 



6 Projection of Pcirticle-Number 



Number-projection is a simple example of the projection theory since in this 

case the matrix Dg is just a phase factor, a multiple of a unit matrix. The 
projection operator which projects out the good particle number [N) is of the 
form 

pN = 1_ /rf0e^<^(^-^). (84) 
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The rotation in the gauge-space is given by 



(85) 



Using this definition in Eqs. (59) and (60), wc obtain the expressions for the 
rotated-densities in the case of particle- number projection as 

= p, e'^t'K. (86) 

This simphfies the expressions for all the matrices defined in section IV. In 
particular, we have 

C^^ ^ pp^ - = e-^V (e^V + - p)) = pC-\<t^) (87) 

with 

C{ct>)^p-'C^^{p-re-'^\l- p))'" (88) 

= e'^'^(l + p(e'^'^-l))~' (89) 

and 

Ct(0) = C(-0). (90) 



For the transition densities we obtain 



P(0) 
k(0) 

«(0) 



:C(0)p, 
:C(0)/€ 



2i<p 



(91) 
(92) 
(93) 



and 



p(-(/)) = Ct(0)p = pt(^)^ (94) 

K{-(f)) = e-2^'^K(0) = C\(I))k = kC*(0), (95) 
K{-(f)) = e-^"t'Ki(f)) = e-'''^«C^((/)) = e-^''^Ci(f))K (96) 



Considering that detD,^ = e we find for the norm-overlap 

x{<P) = i-e-^^^Vdet(e2^V+(l-p)) = ^4=^- ^9^) 

27r ^^let C{(f)) 
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For the Hamiltonian-overlap, we have 



HM^Tr{ep{cf>)), (98) 



2 ^ 2 

711712713714 



HM) = -^Tr (A*((/))«:(0)) = -^IV (A(0)r (0)) , (100) 

(0), (101) 



4 

711712713714 



with 



711713 ^711712713714/^714712 

(0), (102) 

712714 

^nin2i<P) = \l2 '>^mn2n3n4l^n3nA^), (103) 

713714 

^n37^4('^) = ^ 51 <lr^2(<^)^r^ln2n3n4' (104) 

711712 



with 



r(-(/.) = rt(</.)^ (105) 

A(-0)=e-2^'^A(0), (106) 
A(-(/))=e-2^'^A((/)). (107) 

Summarizing, we can write the projected-energy as 

E"" ^ j d(t> y{(t>) (HsM + + HM)) ■ (108) 

In the following subsections, we shall evaluate the variation of the norm- and 
the Hamiltonian-overlaps. 

6. 1 Variation of the Norm 



Using (14), the norm in the case of the particle-number projection can be 
rewritten in terms of the density (p) only. This simplifies the variational equa- 
tions considerably. Defining the matrices X{(f)) and Y{(j)) by 



dx{(t)) 

dpn'n 



(109) 
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^=l/(0)i;.,(0), (110) 

OPn'n 

we find by differentiation 



--a;(0)(l-e-2^<^)C„,,(0) = a;(0)ze-*^ sin C,„,(0), (112) 

and obtain 

X(0) =ie-^'^sin0C((/)), (113) 

and 

F(0) = ^e-*^sin0 C7(0) - ^ y d0V(0')e"'^' sin0' C{(P'). (114) 
5.^ r/ie Projected Hartree-Fock Field 



The HF-field in the projected HFB-equtions is obtained as the derivative of 
the number-projected energy with respect to the density 



C-l^ (115) 
= J d<Py{<P) (Yr,A<P)H{<P) + for (n < n'), (116) 



which gives 



f)F^ r 

h'' = -^ = Jd<P y{<t>)Y{<t>){H,,{<t>) + HM) + HM)) (117) 
+ 1 ci0 y{(t>)— {H,,{cj>) + H,u{4>) + H,^m . (118) 
The above projected HF-potential has three parts 

/,^^£^ + r^ + A^ (119) 



with 
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T'' = jd^ 1/(0) ('F(0)i/,,(0) + ^H,,{<P)\ , 



dp 

A^ = /d0 y(0) (y{cI>)HM) + ■ 



Using 



Tr(A^^C{<P)\^-2ie-''f'sm<l> (C(0)AC(0))„„, , 



we find 

dH,p{<p) ^ f d 



Tr e- [C{(f))p] for (n < n') 



-"n'n 



= {[l-2ie-''f' sin (j)p{(j))]eC 
^^ = Trfr(0)-^[C(0)p]) for (n<n') 

= ([l-2ze-^sin0p(0)]r(0)(:7(0)) 

_le-2i'^TrfA(0)-^k*C]l for (n < n') 

2 V / 

= ie-*^sin0 (/«(0)A*(0)C(0) + (:7(0)A(0)7«*(0)) 
Considering Hermiticity, we finally have 



£^ = ^|#y(0)(y(0)ivM0)] 

+ [1-2^6"*'^ sin 0p(0)]eC(0)^ + h.c. 

r"^ = ^/c^0i/(0)(n0)^^[r(0)p(0)] 

+ i[l-2ie-^'^sin0p(0)]r(0)C(0)') + h.c. 



A^ = -^/# ?/(0)(n</')^T>[A(0)^*(0)] 
-2ie-^*sin0C(0)A(0)7«*(0)) + h.c. 
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6.2.1 The Projected-Pair Field 

The pairing-field is obtained by variation of the projected-energy with respect 

to K 



^nn'--^ for {n<n') (133) 



(134) 



We find 



^^ = -ie-^^Tr(c(0)A(0)^/.*) for [n < n') (135) 

^-\e-'^f(C{<t>)^{<t>)-{..f)^^,, (136) 
which finally yields the pairing-field in the projected HFB-equations 

= i / #e-2^^ y{ct>)C {ct>) A(<l>) - (..f (137) 



7 Canonical representation 



In this section, we shall demonstrate that the projected expression for the 

A-matrix given in the last section reduces to the familiar form in the canon- 
ical basis. In the canonical basis, which we shall denote by the greek indices 
(/x,i/, ...), the density-matrix and the pairing-tensor reduce to the following 
(2x2) matrices 



Using these basic expressions, the matrices in the number-projection, Eqs. 
(89-92) acquire the following form in the canonical basis 



+ (139) 





"1 




( 


Uu,Vu, ^ 


















-Uf,v^ J 
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n K + (140) 



2^ .>o 



p^{ct>)^e^'^CMvl^I (141) 
K^{(t>)^e'^C^{(t>)u^v^xX (142) 
«;(0) = e-''^C^(0)xi^^;^. XT, (143) 



where C, = e^'^'^ and 



1 \ / 1 , 

and X = I . (144) 

1/ I -1 



The norm is therefore given by 



2-K 

< $|P^|$ >= / #a;(0) ^ — UC C""-' n + = K (145) 

where the residuum interals are defined in analogy to ref. [12]: 

K{lx,,...,lXrr.) = ^jdCC"-' n « + C<)- (146) 

/U7^/il../im 



Using the expression for the pair-gap (137) and the canonical forms of the 
matrices in (138-143), we obtain 



YlACul + Cvl) 



iCui + cvDiCui + cvi) 

- (/X ^ fi') (148) 

The diagonal form (^ = /j,') of this expression agrees with the expression 
derived in ref. [12]. The expresions for the other fields in the canonical repre- 
sentation also agree with the expressions given in ref. [12]. 



8 Model study 



As a case study for the projection formalism developed in the present work, 
we use the deformed shell model Hamiltonian which consists of a deformed 
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one-body term, h and a scalar two-body delta-interaction [19]. The one-body 
term is the familiar Nilsson mean-field potential which takes into account of 
the long-range part of the nucleon-nucleon interaction. The residual short- 
range interaction is specified by the delta-interaction. The model Hamiltonian 
employed is given by 

H = -AkJ—Y2o - g5(h - h). (149) 



We use G = g J R^ir^dr as our energy unit and the deformation energy k 
is related to the deformation parameter /3. For the case of hn/2 shell, k,—2A 
approximately corresponds to /3 = 0.25. The model Hamiltonian has been 
solved exactly for particles in /iii/2 intruder subshell and the results will be 
compared with those obtained using the projection formalism. In the numerical 
calculations, we have used only J = and 2 components of the 5-interaction. 
The problem with using the full (5-interaction is that in the HFB analysis of this 
simple model, the contribution of the higher-multipoles to the particle-hole 
channel is quite large and it renormalises the one-body potential substantially. 
This renornialization makes it very difficuilt to choose a reasonable initial 
HFB fields. In most of the mean-field analysis this particle-hole channel is 
not considered. However, since we are comparing the results with the exact 
calculations, it is important to keep all the components of the mean-field. 
Nevertheless, the main purpose here is to choose a simple model which can 
solved exactly and the projection method will be tested using these exact 
results. 

In the HFB analysis of the single-j shell, the single-particle basis are the 
magnetic-states of the single-particle angular- momentum, j. The summation 
indices in all the expressions given in section VI run over these magnetic- 
states. For the case of /111/2 orbital, the summation indices have the range 
(—11/2, ....,11/2) and since the Hamiltonian in (149) obeys the time- reversal 
symmetry, it is required to have only -|-ive or -ive magnetic-states. Therefore, 
the summation indices have dimensions of six for the /iii/2 orbital. 

In order to test the projection method, we have performed HFB, PHFB and 
the exact shell model calculations as a function of the pairing strength, G for 6- 
particles in hii/2- In Fig. 1, the results of total energy (Etot) and pairing energy 
(Epair) are presented. In PHFB we have used six-mesh points for integration 
over the gauge-angle. It is found that for the value of G close to 1, the three 
mesh-points give quite accurate results. However, for the limiting value of G 
close to 0, it is important to use six-mesh points. The PHFB results of Efot 
with six-mesh points reproduce the shell model results almost exactly. In Fig. 
1 the results of PHFB and the shell model calculations are indistinguishable. 
As expected the results of PHFB and HFB are identical for G = since there 
is no two-body interaction and Efot is equal to the energy of the one-body 
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static-potential. The results of PHFB and HFB deviate with increasing value 
of G and the deviation is more than 1 unit for the value of G = 1. The results 
of Epair are shown in the lower pannel of Fig. 1 and it is clear that HFB shows 
a phase transition G — 0.4. For G — 0.4 and lower, -Epair is exactly zero. 
The PHFB on the other hand depicts no such phase transition and has a finite 
value for finite G. 



9 Summciry 

In order to obtain an approximate solution of the many-body problem, it is 
often required to break the symmetries which the original many-body Hamilto- 
nian obeys. For instance, in the case of Hartree-Fock or Hartrcc-Fock-Bogoliubov 
approaches, the many-body wavefunction is approximated by a Slater-determinant 
which in a general case breaks the rotational and the particle-number symme- 
tries. In order to have a better description for the many-body problem, it is 
essential to restore the broken-symmetries. 

The symmetry restoration may not be very critical in strongly symmetry 
breaking conditions, but for the case of a weakly broken symmetry the restora- 
tion is quite important. For example, in the case of well deformed nuclei in the 
ground state, the pair-correlations are strong and the BCS theory which breaks 
the particle-number symmetry, describes the pair-correlation quite reasonably. 
However, with increasing quasiparticle excitation, the pair-correlations reduce 
and the BCS solution collapses. The pair-correlations are also quenched for 
superdeformed shapes. For these weak pairing cases, it is quite essential to con- 
sider the restoration of the particle-number. As already mentioned in the intro- 
duction, the most undesirable feature associated with the mean-field theory is 
that one obtains unphysical phase transitions, for instance from finite-pairing 
to zero-pairing. These phase transitions are smeared out with the restoration 
of the broken symmetries. 

In most cases, the projection methods have been carried out approximately 
[15]. A self-consistent description of the projection in the HFB framework 
has been an unresolved problem. In the present work, we have obtained the 
symmetry-projected HFB equations for the first time. This has been possi- 
ble by first realising that the projcctcd-cncrgy can be completely expressed 
in terms of the HFB densities and then it has been demonstrated that the 
variation of an arbitrary energy functional which is completely expressible in 
terms of HFB densities, results in HFB equations. The expressions for the 
Hartree-Fock and the pair-field depend on the form of the energy functional. 
In the case of the HFB energy, the expressions for these fields are quite sim- 
ple, whereas with the projected-energy the expressions for these quantities are 
more involved. 
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We consider that the major advantage with the present projection formahsm is 
that one can use the existing HFB computer codes and the only the expressions 
for the fields need to be redefined. We have applied the projection method to 
an exactly soluble deformed single-j shell model with the conclusion that the 
numerical effort involed is similar to performing the bare HFB calculations. 
The results of the projection method almost agree completely with the exact 
shell model calculations and the phase transition obtained in HFB is smeared 
out with the projection. 



10 Appendix: Helpful Formulae 

Here, we present some helpful formulae for calculating the derivatives for ma- 
trices and matrix functions which have been extensively used in the present 
work. If A{x) and B{x) are arbitrary matrices depending on a parameter x, 
we have for A' — dA/dx 



d_ 
dx 




A'B + AB\ 



-A-^A'A-^, 



(152) 



(150) 

(151) 



d 

dx 



det(^) 



detylTr(yl-^yl'), 



(153) 




(154) 



dx 




dei{AB) Tr {AT^A' + B-^B') 
dei{A-^B) Tr (-A'^A' + B'^B') 



(155) 



(156) 
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Fig. 1. The results of the total energy (Etot) and the pair energy (Epair) for 
six-particles in a deformed j = 11/2 orbitial using different approaches. The pro- 
jected-HFB results for Etot are indistinguishable from the exact shell model calcu- 
lations. 
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